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FOREWORD 


This report summarizes a three-year scope of work at the University of 
Akron within the field of transient/nonlinear dynamic analysis of gas 
tnrhine jet engine strnctnres. Dr. C. C. Cham is was the NASA-Lewis 
Technical monitor of this grant. The focns of this work was on the 
development and implementation of rotor-stator interactive force elements 
for implant into general purpose nonlinear time-transient finite-element • 
codes snitable for general engine dynamic simulation. The Principal 
Investigators on this grant were Drs. K, L. Adams*. J. Padovan and D. G. 
Fertis of the University of Akron. Graduate students involved with various 
portions of this work were I. Zeid, R. Quinn. E. Eirt. S. Payandeh and 
P. Lam. 


* Now Professor of Mechanical Engineering at Case Western Reserve 
University. Cleveland. Ohio 44106. 



1. INTRODUCTION AND BACKGROUND 

There is currently a considerable interest and level of activity in 
developing computational schemes to predict general engine dynamic 
behavior. The general feeling among knowledgeable parctitioners working on 
engine vibration problems is that various modes of operation, such as 
blade-loss events, require a high level of analysis sophistication to 
realistically model the engine. Proper account of system nonl inearities 
(particularly at the bearings, dampers and rubs) appears to be necessary if 
analytical predictions are to be realistic. The approach used in this work 
makes use of already proven general finite-element nonlinear time-transient 
computer codes which are available on the optn market. The work specifi- 
cally described covers a three-phase NASA-Lewis sponsored research grant on 
engine dynamic simulation with available finite element codes. Ptm — A 
concentrated on the development of a bearing-damper element computer soft- 
ware package suitable for 'plug-in* to available finite element codes. 

Ph,te -2 focused on the implanting of this damper element into the ADINA 
code and extensive computational testing of the completed package. EfetSi zl 
was related experimental work on a new test rig employing a flexible rotor 

nountod in two sQDCCZC"filB dampers* 

Present-day jet engine configurations have evolved to a substantial 
degree through a trial-and-error process involving extensive testing. There 
are many fundamental dynamic phenomena which take place within these engines 
for which basic description and understanding have yet to be generated. 
Nonetheless, they work well. Modern aircraft engines are typical of current 
high-technology products in which the recently acquired computing capabili- 
ties of today are being used to better understand and improve what is 
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already designed, built and operating. 

A better understanding of tbe basic dynamic characteristics of existing 
and nev engine configurations is a prerequisite for producing acceptable 
engine efficiencies on advanced configurations (i.e., smaller rotor/stator 
running clearances). Also, a better definition of engine dynamic response 
would more than likely provide valuable information and insights leading to 
reduced maintenance and overhaul costs on existing configurations. Further- 
more, application of advanced engine dynamic simulation methods could 
potentially provide a considerable cost reduction in the development of new 
engine configurations by eliminating some of the trial-and-error process 
done with engine hardware development. 

The emergence of advanced finite element codes, such as NASTRAN, 

NONSAP, MARC, ADINA, ANSYS and ABAQUS and related algorithmic advances, has 
placed comprehensive engine system dynamic analyses within reasonable reach. 
Vhat remains to be done is to develop new component element software to 
properly model engine rotor/stator interactive components, such as the 
squeeze-film damper, within tbe algorithmic logic of already proven finite 
element codes. This has been a major mission of this work. 

For good reasons, aircraft gas turbine engines use rolling element 
bearings exclusively. This design philosophy had, in earlier times, 
deprived engines of the beneficial damping inherent in many other types of 
rotating machinery where fluid-film journal bearings are used. The 
implementation of squeeze-film dampers in later engine designs has now 
provided engine designers with an effective means of vibration energy dis- 
sipation. The net result is that engines with squeeze-film dampers are less 
sensitive to residual rotor imbalance and better able to control vibration 
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and transmitted force levels resulting from various excitation sources 
within the engine. 

The field of rotor dynamics has evolved to its present state primarily 
through the solution to problems in classes of machinery older than aircraft 
engines. In most other types of rotating machinery (e.g., steam turbines, 
centrifugal pumps and compressors, fans, generators, motors, etc.) the rotor 
can be adequately modelled as an Euler or Timoshenko beam [1]. In addition, 
the support structure holding each bearing can often be adequately modelled 
as a separate mass-damping-stiffness path to ground (i.e., to the inertial 
frame). Also, for most purposes, bearing dynamic properties are 
characterized as stiffness and damping elements, linearized for small 
vibration amplitudes about some static equilibrium state. Tith few excep- 
tions (e.g., Hibner [2]), it is this level of sophistication that has been 
utilized for the most part is rotor-dynamics analyses of aircraft engines. 

Present day aircraft engines are structurally far more complex than 
most other types of rotating machinery. The multi-shaft configuration, plus 
the fact that the shafts are thin rotating shells, creates unique but sig- 
nificant complicating differences between aircraft engines and otber turbo- 
machinery. Also, the stator structural support at each rotor bearing repre- 
sents anything but a separate mass-damper-stiffness path to an inertial 
frame. In fact, setting the inertial frame for the engine is not a simple 
matter when the full range of in-service maneuvers is realized. Dynamic 
paths between different bearings exist not only throogh the rotor but 
through several other paths within the nonrotating engine structure, i.e., a 
'multi-level', 'multi-branch' system. As many as eight significant 'levels' 
have been identified. 
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The feasibility of nonlinear dynamic analyses of multi-bearing flexible 
rotors has been recently demonstrated on non-aircraft applications 131. 
There are highly nonlinear dynamic effect in aircraft engines, particularly 
under large excitation forces, such as blade or disk failures, hard land- 
ings, and foreign matter ingestion events. 

Clearly, the field of aircraft engine dynamics is presently in a 
position where there is both a need for substantial advances and feasible 
means available by which such advances can be accomplished* 
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2. TIME-TRANSIENT NONLINEAR DYNAMIC ANALYSES 


In recent years it has become evident that an important class of engine 
dynamic phenomena can not be studied without accounting for the highly non- 
linear forces produced at bearings/dampers, labyrinths and other close- 
running rotor/stator clearances under large amplitude vibrations. In such 
cases, linear theory typically predicts vibration amplitudes larger than the 
actual running clearances. Furthermore, important vibratory phenomena, such 
as subharmonic resonance and motion limit cycles, are 'filtered' out of the 
problem with a linear model, giving grossly erroneous predictions, 
qualitatively as veil as quantitatively [3]. 

With few exceptions, nonlinear dynamics problems must be solved 
numerically as time-transient responses, whether the sought answer is a 
steady-state periodic motion or is strictly a transient phenomenon. The 
problem is mathematically categorized as an initial value problem in which 
the displacements and velocities of the complete system must all be 
specified at the beginning of the transient. From that point forward in 
time, the equations of motion are numerically integrated (known as 
'marching') as far in time as one wishes to study the system motions and 
forces. If the system is dynamically stable, the transient motion dies out 
yielding the steady state response which in most real systems with a perio- 
dic force excitation will be a periodic motion. In a stable system with no 
time-varying force excitation, the transient will die out as the system 
comes to rest at one of its stable static equilibrium positions. If the 
system is unstable, the transient does not die out but continues to grow in 
time unless or until some nonlinear mechanism in the system limits the 
motion to what is frequently called a 'limit cycle' [3]. 
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In order to study tbe general dynamical characteristics of aircraft 
engines, nonlinear dynamics computational schemes are required* The 
approach taken is to develop software packages to model engine components 
which are not typically found on dynamical structures and therefore are not 
already built into existing nonlinear finite element structural dynamics 
computer codes. The initial effort has concentrated on developing such a 
software package for squeeze-film bearing dampers. 
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3. oyFPAT.T. APPROACH USING INTERACTIVE ELEMENTS 

Considering the typical engine structural complexities, an improved 
computational approach is necessary if a proper transient/steady-state model 
is to be developed for gas turbine engines. In this approach, it appears 
that the finite element method is one of the attractive modelling techniques 
for such problems. Its inherent capabilities include features essential to 
modern engines: 1) automatically handles multi-branch, multi-level struc- 

tures in a more direct and efficient manner than flexibility approaches. 

2) well-suited to handle nonlinearities associated with structural kinematic 
and kinetic effects [4], 3) easily accommodates various types of boundary 
and constraint conditions, and 4) easily accommodates material nonisotropy 
and nonlinearity 14.51. A body of established and proven algorithms are 
available which can handle these various important effects [4.6] as well as 
geometric complexities, e.g., beam, plate, 2-D and 3-D elements [71. 

The required features which were not previously available with general 
purpose finite-element codes are provisions to handle rotor/stator 
interactive forces originating from squeeze-film dampers, seals and rub/ 
impact events. Presented in the next three Sections of this report are the 
results of an effort to develop a squeeze-film damper computer software 
package which can be 'plugged' into any existing finite element code. 
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4. SQUEEZE-FILM DAMPER ELEMENT DEVELOPMEm - 


The bearing damper element is essentially an interactive element to re- 
present squeeze-film dampers. Its purpose is to bridge the 'gap' between 
structural elements which are separated in the actual engine by a squeeze- 
film damper. In its simplest version, it has an input/output setup as shown 
in Fig. 1. A source listing of this code is given in reference [8]. 

The rotor/stator interactive force generated in a bearing squeeze-film 
damper is modelled using an adaptation of the classical Reynolds lubrication 
equation for incompressible laminar isoviscous films. 


where. 


3 / h' 


9 i 

(£ dp \ 

rjz l 

^ ft dz ' 


dh 


(1) 


z * axial coordinate 
z * circumferential coordinate « rO 
h * local film thickness 
dh/dt ■ instantaneous local rate of change in b 

U * sliding velocity * Rid, typically zero in a damper 
C « radial clearance of damper annulus. 

The relationship between system inertial coordinates and damper parameters 
comes through the expressions for h, dh/9x and dh/dt. Referring to Fig. 2, 
these relationships are summarized as follows: 




( 2 ) 


^(X M -X S )U(Y„-Y,)j 


(3) 


then 


h = C- e-n = C- ( X* - X s ) con 6- ( V* - >\ ) sin 9 


(4) 


8 




dh I dh I 


( 5 ) 


and 

dh ... 

— »-(*,-A\>cosd-(F„-y,)«nd (6) 

Many configuration* of dampers which are currently being employed in 
engines have done away with centering springs common in older designs. 

Design simplicity as well as centering-spring fatigue life are apparently 
the major reasons. Also, in the majority of cases, damper end seals are 
used because this keeps damper throughflow sufficiently low to be compatible 
with the overall engine lub system of pre— damper configurations. However, 
the disadvantage of having end seals is that a 'large' damper clearance of 
typically 10 mils is required in order for the squeeze-film action to effec- 
tively dissipate vibration energy. From other considerations, a smaller 
damper clearance would be desirable (e.g., blade tip clearances). Without 
end seals, the optimum damper clearance is considerably smaller. Engines 
with higher oil flow capacity and no damper end seals are probably the trend 
on future designs. 

Two typical configurations are shown in Figs. 3 and 4. The end seal 
configuration in Figs. 3 and 4(a) essentially divides the lubricant annulus 
into two pressure domains whereas that in Fig. 4(b) is a one-domain problem. 
In both cases, the 'long-bearing' solution is appropriate. Both solutions 
are options in the software package developed in this work. 

For the 'long— bearing' solution dp/dz (( dp/dx, and the following ordi- 
nary differential equation (two point boundary value problem) is obtained 
from equation (1), for U * 0. 
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For configurations with no end seals, an improved adaptation of the 
short-bearing approach is used by implementing the parabolic assumption of 
O'Donoghue [9]. The following approximation is made: 


P(0.C)=p(0,O)(l- 



( 8 ) 


This assumes an axially symmetric axial pressure distribution at every cir- 
cumferential location, and results in the following pressure field equation. 



o piB % 0)h' 


(9) 


This is actually a first-order Fourier approximation using the parabola as 
the single approximative function. 

A convergent approximation to the full two-dimensional Reynolds equation 
can be obtained, as an extension of the foregoing approach by O'Donoghue [9]. 
The number of Fourier terms is increased to N, resulting in N simultaneous 
ordinary differential equations. 


p ( 6.Z )*=/>, (0.0) cos ^ + P: (0.0) cos ^ + . 


*P\ (0.0) cos 




( 10 ) 


Substitution into the general 2-D Reynolds equation (1), followed by LHS:RHS 
segregation by arguments yields N ordinary differential equations, one for 
each p. (0,0) , [9] . 


5. METHOD fiE DAMPER ELEMENT SOLUTION 


Although there are computationally fast closed-form solutions available 
such as given in [10,11], they do not retain sufficient generality to handle 
the specified pressure boundaries at supply and drain ports (Figs. 3,4,5) of 
typical configurations. Also, they are not amenable to structural deflec- 
tions of the damper elements. For these reasons, the following approach has 
been nsed. Although somewhat computationally slower than closed-form solu- 
tions, it is nonetheless computationally quite efficient, and retains the 
generality deemed necessary. The solution method given below is used for all 
three formulations, i.e., equations (7) through (10). It is described below 
as implemented for the long-bearing formulation. 

Based on a 3-point central difference, the following long-bearing equa- 
tion yields equation (11). 


d\ d\ dx dr 


< i ). 


m p. . i - p , , 


dx > ■ 2Ax 

d-p \ _ P,. t -2P,+ P, , 




Av : 


K av / dx v & 


Rearranging (11) produces 


( 11 ) 



_ 3A? dM dh, 

. wl 1 i 2±x dx \ ** dt 


which is condensed to the following form: 

C,P, + tjP, , ■* D.P,., 


( 12 ) 
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Employing the recursion relationship 


P ,- 1 = A,P,+B, 


( 13 ) 


equation (12) can be expressed in terms of only two adjacent grid points as 
follows : 


or 


Therefore , 




(14) 


P J {C,+E J A l )+E 1 B t + D,P l ^=R, 


(15) 


p ~ ( Bi ^/> + zJisBi 

' \C t -iE,A t ) *' C,+E,A, 


(16) 


which, when compared to equation (13), yields the following recursion 
relationships: 


A,, i * - 


C,+E,A, 


B - R ~ EB ' 


(17) 


(18) 


c,+e,a, 

From the upstream boundary condition for each domain, the {A} and (B) vectors 
are determined by starting with Aj * 0, Bj * P^ (called the forward sweep). 

The downstream boundary condition is inserted at the beginning of the 
backward sweep as follows: 


Pm- i m A mP \t + B\t 

i+ Aw i 

P; A \P 1 B\ 


(19) 


Film rupture is handled by the following substitution. If Pj < P va p 0 r* * et 
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Pj * Pyapor ^ c ^ orc computing Pj.j, This is equivalent to the condition dp/dx 
* 0 at the film-rupture full-film boundary. In the case of the 2-D 
convergent approach indicated by equation (10), this point-by-point test is 
made on the local summation. 

PH'**- r pa9, . 0 . 

A I 

The method of solution, although not closed-form, is noniterative and is 
a 1-D adaptation of the 2-D finite difference method of Castelli and Shapiro 
[12]. While it does entail a one-dimensional, finite-difference scheme, it 
requires only a very small amount of CPU time and is therefore ideally suited 
to time transient rotor dynamics analyses. It has major advantages over the 
purely closed-form approximations, e.g., [10,11], as noted earlier. These 
major advantages are immediate account of specified-pressure boundary condi- 
tions at feed and drain holes of a damper. Also, the finite difference 
approach easily permits account of static as veil as dynamic deflections 
which alter the oil film gap geometry from ideal rigid circular shapes. 

Forces components on rotor are computed by numerical integration of the 
instantaneous film pressure distribution, as is standard. 


= - \ cos SdA=-LR\ 'pWcosMB 

' J 4f J "| 

F, = - j r i\nBdA= -LR^'p{e)i\n»de 


( 20 ) 


Stator force conponents are equal but opposite the rotor force 
components. 

Force gradients (i.e., instantaneous tangent stiffness and damping) com 
ponents are obtained by local 'small' perturbations, as is standard. 





( 21 ) 


IS 


Numerical differentiation is performed with small AXj and AXj increments 
about instantaneous conditions. This provides continuous updating of {Fj} f 


6. pEMON STRAT ION 0£ DAMPER ELEMENT 

For purposes of checking out the damper element code (SQUEEZ) and to 
demonstrate its use, two types of computations were made. First, a paramet- 
ric study of damper pressure distributions was made for a variety of speci- 
fied circular orbits, for both long-bearing and short-bearing solutions* 
Second, a four-degree-of-freedom rotor-damper-stator model was investigated 
under conditions of small rotor unbalance through large rotor unbalance* 

These results are summarized* 

For the parametric study on pressure distribution, the following damper 
annulus parameters were used: 

Diameter, D = 6 in. 

Length, L * 1*25 in* 

Radial clearance, C * 0.010 in. 

Lubricant viscosity. |i « 1 * 10" 6 reyns 

Angle between inlet oil port and drain port, (Oj - 0 O ) « 180 deg 

Inlet oil port pressure, Pj * 55 psia 

Drain port pressure, p c * 15 psia 

Lubricant vapor pressure, p v * 1»5 psia 

Orbit angular velocity, 0 - 3600 cpm (377.rad/sec) 

The above damper parameters are typical for modern gas turbine aircraft 
engines. A parametric stndy was made postulating the outer ring of the 
damper fixed and the inner ring having a constant-radius constant-velocity 
concentric orbit. Eccentricity ratios (i.e., orbit radins/radia 1 clearance) 
from 0.05 to 0.95 were computed, both for the long-bearing and short-bearing 
solntions. 

Circumferential center-line pressures were plotted as a function of cir- 
cumferential position and time, for one period of prescribed motion. The 
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results for the long-hearing solution are shown in Fig. 6, and for the short- 
hearing solution in Fig. 7. The difference between long-hearing and short- 
hearing solution is quite large when compared with the same radial clearance. 
One therefore sees why dampers with end seals require larger clearances to 
work properly than dampers without end seals. 

A simple 'driver' code was written (see listing [8]) which uses the 
damper-element code in the same manner as a general implanted application 
with large finite element codes. The 'driver' code is based on a four- 
degree-of— freedom system, i.e., planar motion of the inner and outer damper 
elements. This then simulates a single-mass rotor connected to a single- 
mass stator through the damper element. The system analyzed is shown in 
Fig. 8. The model is coded to simulate arbitrary rotating and/or static 
radial loads. Aside from demonstration purposes, this four-degree-of- 
freedom model has been devised as a 'bench tester' of different damper 
models and configurations. 

Note from Fig. 8 that the high-pressure port (i.e., feed port) is 
located on the bottom of the damper so as to assist 'lift-off'. Since 
centering springs are not typically used, they have been excluded in this 
example, being the most nonlinear type case and therefore computationally 
the most demanding of the algorithms employed. Lift-off therefore requires 
some amount of vibration to overcome the dead-weight load. Rotating 
unbalance loads of 100, 200, 300, 500, and 1000 lbs were run with 0 • 150 
rad/s. Orbital plots were made showing rotor and stator total motion on one 
plot and rotor-relative-to-stator motion on a second plot. The plotted 
results are shown in Figs. 9 through 13. 

For a 100 lb rotating load (Fig. 9) the motions shown are for a 20 
load-cycle transient from time * 0. The rotor and stator each show close to 



the same motion, end their relative motion is small, with the rotor barely 
'lifting off'. The relative orbit is essentially oscillatory. However, 
when the rotating load is increased to 200 lbs, (Fig. 10), the relative 
orbital motion shows the beginnings of orbital motion, i.e., a 'crescent 
moon' shape as measured by numerous investigators. Further increase in 
magnitude of the rotating load to 300 lbs (Fig. 11) shows a well-defined 
steady-state total motion as well as relative motion. Note that with a 300 
lbs rotating load, the relative (rotor-to-stator) orbit is still small in 
comparison to the radial damper clearance and confined to the region of the 
bottom of the damper. However, an increase of rotating load magnitude to 
500 lbs causes a considerable change to the relative orbit (Fig. 12). 

Notice now that the relative motion of the rotor with respect to the stator 
fills a major portion of the clearance circle. Further increase of rotating 
load magnitude to 1000 lbs (Fig. 13) simply causes the steady— state relative 
orbit to expand and fill even more of the damper clearance circle. 
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7. SQUEEZE-FILM DAMPER IMPLANT INTO ADINA - A DEMONSTRATION 


The squeeze-film damper element (SQUEEZ), described in previous 
Sections of this report, is also summarized in references [8] and [13], 
Reference [8] also contains Fortran computer code listings for the damper 
element as well as for the four-degree-of-freedom driver code used in its 
development and demonstration. 

To demonstrate the implementation of the SQUEEZ element into a general 
purpose transient nonlinear finite element (FE) code the ADINA code was 
purchased from MIT by the University of Akron as a cost-sharing expenditure. 
The next group of Sections in this report describe the implant of SQUEEZ 
into the ADINA code, a general purpose nonlinear FE code well known and used 
throughout the field of structural mechanics. This implant work is also 
summarized in references [14] and [15], with complete Fortran listings of 
all generated software and graphics packages given in reference [15]. 

To benchmark the approach described in references [14,15], the appro- 
priate general purpose finite element code had to be selected. As noted in 
the first year report [13], the ADINA code was chosen for this purpose. 

This follows from the fact it has the requisite features required for rotor- 
bearing-stator simulations. These include: (i) Linear and nonlinear 
substructuring features, (ii) An extensive element library, (iii) Capability 
to handle kinetic, kinematic and material nonlinearity, (iv) Explicit and 
implicit integration loops, (v) Simplified I/O features, (vi) Accessible 
code architecture, (vii) Extensively benchmarked, and (viii) Requisite 
portability and general availability. In terms of the squeeze-film damper 
element noted in [8], extensive modifications were introduced into the ADINA 
code. These modifications were made general enough so as to handle rotor- 
bearing-stator simulations involving any number of rotors and associated 
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squeeze-film damper elements. The ADINA program has two available solution 
procedural loops, namely either explicit or implicit time-step integration 
of the equations of motion. 

Since extensive amounts of data are generated during a typical run, 
graphics post processors have also been developed to simplify output eva- 
luation. These include both 2-D as well as 3-D plotter schemes. The 2-D 
processor enables the plotting of: (i) Bearing/rotor trajectories at a 

given station, (ii) Stator trajectories, (iii) Clearance histories at 
given bearing stations, and (iv) Force, velocity and acceleration histories 
at given bearing stations. The 3-D processor enables the plotting of 
isometric views of the rotor trajectories. All these code listings are con- 
tained in reference [15]. 
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8, NUMERICAL CONSIDERATIONS AND EXAMPLES 


This Section describes the work pertaining to the computational charac- 
teristics of ADINA-with-SQUEEZ direct integration approach of FE generated 
rotor-bearing-housing motion simulations, 

8.1 Benchmarking 

To benchmark the overall procedure, a simple lumped parameter direct 
integration scheme was developed. This approach was used to check the 
accuracy of the FE generated scheme involving the ADINA-SGUEEZ 'implant*. 

As the first example of such benchmarking, consider the system defined in 
Figure (14). The material and geometric properties associated with the dual 
ported squeeze-film damper bearing employed in this and the following sample 
problems are defined by 

Nominal diameter *= 6 inches, 
nominal length * 1.2 inches. 

Clearance * 1. inch, viscosity * .1*10^, 
film rupture pressure * 15 psia, 

* 90, ©2 * 270, Pj * 15 psi, Pj * 55 psi. Note as with current 
practice, the bearing used in the simulation has no centering spring. In 
terms of this system. Figs. (15-17) illustrate various aspects of the 
validation of the Newark [16], Vilson [17] and central difference generated 
results with those of the benchmark. In this way, both the implicit and ex- 
plicit schemes are treated. In each of the comparisons, three different 
aspects of rotor bearing stator behavior are depicted specifically. 

a) Rotor displacement trajectories, 

b) Stator displacement trajectories, and 

c) Relative rotor orbit. 
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For the present purposes, while extensive benchmarking was undertaken, 
for convenience, only the case of mild rotor unbalance is depicted. As can 
be seen, for the given At steps chosen, good comparisons were obtained by 
both the implicit and explicit schemes. Such benchmarking was obtained over 
a wide range of rotor speeds and rotor unbalance levels. In all cases, good 
accuracy was yielded. Similar benchmarking was also performed for multi- 
bearing-rotor-stator simulations. Note, so long as At was kept small, good 
accuracy was obtained over a wide range of system variables. 

8.2 Explicit and Implicit Direct Integration Methods 

For problems involving few degrees of freedom, it was found that for a 
given accuracy, both the implicit and explicit schemes required about the 
same overall computational times. This follows from the fact that for 
’small problems', the architectural overhead associated with the implant 
strategy programming dominates over the relative algorithmic efficiency. 

Note for transients initiated by rather severe unbalance loads, it was found 
that the implicit scheme proved to be more sensitive to the choice of time 
step size. Interestingly, such, sensitivities were found to occur for 
problems with small as well as large numbers of degrees— of— freedom. After 
performing several parametric studies, it was found that during the course 
of a typical transient, particularly involving a severe loading, the tangent 
properties of the fluid film undergo major changes. Because of this, a per- 
fusion of system harmonics are introduced into the transient response as the 
tangent fluid film properties vary. Note, this behavior is intrinsic to the 
fluid film and hence, is independent of the number of degrees-of-freedom of 
the rotor-stator model. Such properties tend to reduce the stability thres- 
hold of the implicit scheme which is best employed when only a few harmonics 
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arc excited. For problems involving large numbers of degrees-of-f reedom, 
the sensitivity of the implicit scheme coupled with the required continuous 
updating and inversions of the dynamic stiffness tends to reduce the running 
efficiency of the procedure. In contrast, the central difference approach 
tends to be less sensitive to such tangent property fluctuations. This is 
true so long as the resulting family of harmonics is bound by the spectral 
characteristics of the rotor— stator system. Because of this, similar At can 
be employed by the implicit and explicit scheme. In this context, the 
explicit scheme has a somewhat better computational efficiency than the 
implicit approach. This follows from the fact that no continuous inversion 
is required by the explicit scheme. 

8.3 Demonstration Problems 

To demonstrate the various capacities of the ADINA implant strategy, 
the results of several example problems will be considered. This will 
involve single and multiple bearing problems with various types of rotor 
unbalance histories, impact events and rotor speeds. For example, in terms 
of the single bearing system given in Fig. (14), Figs. (18 and 19) 
illustrate various aspects of the response histories to an imbalance load 
which is applied as a ramp function in time. As this loading is more severe 
than that applied in Figs. (15-17), the rotor tends to fill its clearance 
circle as seen in Fig. (19). If the same system is subject to a unidirec- 
tional impulse, as might be expected from a rough landing, Figs. (18) and 
(20) illustrate the associated response history. By comparing Figs. (18) 
and (20), the effects of the impulse can be clearly seen from the ovalixing 
of the rotor trajectories. This is directly due to the directional charac- 
teristics of the impulse load. Note, comparing the results in Fig. (19), we 
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see that the rotor orbits are essentially unchanged by the presence of the 
shock load. Hence, it follows that the squeeze-film damper has essentially 
no effect on mitigating the worst aspects of nnidirectional shocks when 
large unbalances are simultaneously present. 

The next series of examples pertains to the rotor-bearing-stator system 
defined in Fig. (21). To simulate the shafting, the beam elements available 
in the ADINA system are utilized. Note, the mass effects are handled via 
the lumped parameter approach. Figs. (22) and (23) illustrate the effects 
of increasing the severity of loading on the dynamic response of multi- 
bearing problems. As can be seen, increasingly stiffer squeeze-film damper 
responses are excited. Note, as the load is further increased, the rotor 
stator tranjectories become 'locked in'. The next series of Figures, namely 
(24-25), illustrate the effects of suddenly applied rotor unbalances. Since 
the rotor speed considered is high, only a small portion of the clearance 
circle is filled. This follows from the fact that due to the rather rapid 
changes in orientation of the exciting load, inertia filtering occurs. 
Because of this, the rotor is supported by severe velocity gradients which 
are only in a close neighborhood. Hence, the pressure gradient generated by 
the inlet and outlet ports of the squeeze-film device causes the rotor to 
settle in the direction of the low pressure port. Note, as the rotor speed 
is decreased, increasingly larger rotor orbits occur. Similar trends occur 
as the level of unbalance is increased. 

Based on numerous parametric studies involving systems similar to the 
foregoing, it was found that all the time integration schemes considered 
were stable for situations wherein the fluid underwent only moderate changes 
in stiffnessduring the overall cycle. Eventhough a perfusion of harmonics 
is introduced by even moderate changes in stiffness (so long as the result- 
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iog spectra are strongly bound by the frequency envelope of the dominate 
system frequencies), spurious energy flow to higher order modes is insigni- 
ficant. Specifically, for the implicit scheme, if the choice of time step 
size is gauged to the dominate higher order system frequencies, then the 
introduction of lower order spectra by the squeeze-film has little effect on 
numerical stability. In contrast, if strong stiffness modulations occur, 
the significant amounts of energy flow are introduced in the ever shifting 
higher order modes. This leads to solution instabilities unless smaller At 
are introduced. In this context, the use of the explicit scheme is 
advocated over the implicit for problems involving strong-to-weak unbalance 
loads . 
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9. EXPERIMENTAL STUDIES 


This Section of the report summarizes the third (final) year work on 
this grant. This effort was focused primarily on basic experimental studies 
on the orbital vibration characteristics of nncentral ized (i.e., no 
centering spring) squeeze-film dampers supporting a flexible rotor. The 
objectives of this work were (i) to conceive, design and build a new test 
apparatus involving two squeeze-film dampers supporting a significantly 
flexible rotor, (ii) to perform a wide variety of experimental parametric 
studies with the apparatus, and (iii) make fundamental comparison with com- 
putational results based on mathematical models and computer algorithms 
developed in the previous two phases of this grant. These objectives are 
snnmarized as follows: 

• Conceive, design and construct test rig 

• Experimental parametric study 

• Comparisons with computational results 

The squeeze-film damper (SFD) work summarized in this Section is 
presented with considerably more detail in reference [31]. the MSME thesis 
of R.D. Quinn. 

9.1 Review si Some Recent ££D SlMl&l 

In 1970, Gunter published an investigation concerning the determination 
of the desirable values of rolling element bearing support stiffness and 
damping for rigid rotors [18]. The mathematical model simulated a general 
fonr degree-of-freedom unbalanced rotor mounted on two damped, linearly 
flexible supports. The stiffness and damping coefficients could be constant 
or speed dependent. 
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In 1973, Giberson [19J, of Turbo Research Inc., published a paper ac- 
claiming the ability of SFD to stabilize otherwise unstable rotors supported 
with journal bearings. Vhen rotor speed surpasses the threshold of instabi- 
lity for the journal bearings, 'whirl* amplitude theoretically grows without 
bound, possibly destroying the machine. Pivoted pad journal bearings can 
stabilize a system by themselves, but 'pivoted pad bearings supported by SFD 
is the arrangement with the greatest ability to stabilize a rotor bearing 
system that uses hydrodynamic bearings*. 

In 1974, Hohan and Hahn [20] published a paper on the 'Design of SFD 
Supports for Rigid Rotors'. The publication includes a general design guide 
for centrally pre loaded SFD, supporting a rigid rotor, mounted on rolling 
contact bearings. The authors also investigated the transient solutions of 
the fully non-linear equations of motion to discover the effect of central 
preloading . 

In 1975, Cunningham, Fleming, and Gunter [21] published an analytical 
investigation on the design of a centrally preloaded SFD for a multimass 
flexible rotor mounted on rolling contact bearings. The paper demonstrates 
a technique for the use of single mass, symmetric, flexible rotor analysis 
to optimize the stiffness and damping for a symmetric five mass rotor. The 
single mass analysis was taken from a 1972 NASA Report by Kirk and Gunter. 

In 1979, Tonneson [22], of the Technical University of Denmark, pub- 
lished a study on experimental squeeze-film bearing orbits. The purpose of 
the investigation was to compare experimentally measured damping coeffi- 
cients with those obtained from simple linearized theory. Comparisons were 
made to both eccentric (offset) and concentric (centralized) damper journal 
orbits. 
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In 1977, Bansal and Hibner [23], of Pratt and Whitney, published an 
'Experimental and Analytical investigation of (dynamic) SFD Forces Induced 
by Offset Circular Whirl Orbits'. In agreement with Tonneson, the authors 
suggest that offset orbits are more realistic for flexible damped supports. 

In 1977, Hibner, Bansal and Bnono [24], of Pratt and Whitney, published 
a paper on the control of instability for the intershaft SFD. An earlier 
analytical and experimental investigation had shown the existence of an 
intershaft viscous damper instability. 

In 1979, Hahn [25] published an analytical investigation on the 'Unbal- 
ance Behaviour of Squeeze Film Supported Rigid Rotors'. Equilibrium load 
capacity and transmissibil ity data for a wide range of operating conditions 
are presented for centrally preloaded and nnloaded vertical rigid rotors 
Donntod on rolling element bearings. 

In 1979, Vance [26], of Gas Turbine Laboratories, published a review of 
rotor dynamics which included the current role of SFD as a solution to 
dynamical problems. 

In 1979, Hibner and Bansal [27] published an experimental and 
theoretical investigation into 'The Effects of Fluid Compressibility on 
Viscous Dampers'. The investigation was begun because of the questionable 
validity of the Reynolds equation when cavitation exists in the squeeze 
film. The effects of cavitation could explain why correlations between 
theory and experiment for SFD have been excellent in some cases and poor in 

others. 

In 1979, Holmes [28], of the University of Sussex, U.I., published a 
study on the control of rotor vibration using SFD with and without a 
parallel flexible element. The experimental and analytical investiga- 
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tion includes a general design philosophy. 

In 1980 and 1981, Cookson and Kossa [29,30] published the results of 
their experimental and analytical investigations into the effectiveness of 
SFD, without parallel flexible supports, used for supporting flexible 
rotors. 

9.2 Experiaental Facility 

The rotor system vas designed for use on a Bently-Nevada rotor dynamics 
test rig. The rotor was powered by a 1/10 horsepower infinitely controll- 
able drive motor capable of driving the test rotor at speeds in excess of 
10,000 RPM. The system included a shaft supported by two sets of preloaded 
duplex ball bearings, with each duplex set mounted in an uncentralized SFD. 
A quill shaft vas used to couple the motor and rotor. Rotor discs of 
various weights could be located at any position between the bearing 
stations (see Figs. 26 and 27). 

The preloaded duplex ball bearings provided for stiff bearing-rotor 
coupling, so that the bearings could be assumed to be rigid in comparison 
with the rotor and damper film. The preload on the bearings also increased 
the rolling friction which decreased the maximum attainable rotor speed and 
increased heat generation. 

The ball bearing housings which preloaded the duplex bearings acted as 
the damper journals. The damper journal diameters were 3.3123 inches, the 
journal land lengths were .475 inches. So, the length to diameter ratio 
was .143, suitably short for a short bearing approximation of the Reynolds 
equation to apply. 

The damper housings had interchangeable inserts allowing the radial 
clearance for the fluid annulus between the damper bearings and journal to 
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be varied from .004 to 0.20 inches. Sealing between the jonrnal and damper 
bearing was accomplished with 0-rings in shear. The lateral clearance 
between the damper jonrnal and housing was dimensioned to yield .007 inches 
of 'squeeze' on each 0-ring. This amount of 0— ring squeeze was determined 
with a trial and error approach to be near the optimum, sealing adequately 
under the required oil pressure yet contributing minimal impedance to 
journal motion. Circumferential grooves with outlet ports were located in 
the sides of the damper housings at the edges of the fluid film annulus. 
Ports were also provided to check the static oil pressure in the grooves at 
the end of the annulus with a pressure gauge. Also, inlet and outlet ports 
were located centrally in the damper inserts and housings at the top and 
bottom of the dampers. Flow metering valves at the outlet ports allowed for 
the end seal or no end seal conditions (or any amount of end leakage or end 
pressure desired). On/off valves at the inlet ports allowed for 
circumferential groove or port oil feed into the damper annulus. With the 
appropriate ports closed or open, the damper conf iguraiton could be model- 
led validly by either the short or long bearing approximation of the 
Reynolds equation. An antirotation pin was threaded through the outside of 
the damper housing, through the inserts and into an oversized hole in the 
journal, allowing journal translation but not rotation (see Fig. 28). Also, 
see Figs. 28 through 32 for more details on the various specially designed 
damper and other parts. 

The steel rotor shaft was 13/16 inches in diameter and 12 inches long 
between bearings. The quill shaft used as the flexible coupling between the 
motor and the rotor was 6 inches long and an 1/8 inch in diameter. See 
reference [16] for the rotor and quill sizing calculations. The flexible 
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coupling eliminated the need for exact alignment of the motor and the two 
rotor bearings. Tiro rotor discs could be located at any position between 
the bearings with set screws. Each rotor disc had threaded holes through 
its thickness located symmetrically about its outside circumference for the 
positioning of balancing weights. The discs were modelled after those that 
were furnished with the Bent ly-Nevada test rig. With both discs axially 
centered on the shaft* the first flexible mode could be excited. 

A gear pump was used to pump the oil from a heated bath through the 
SFP. A metered bypass line provided a wide range of damper inlet pressure 
and flow rate. Heating the SAE 10 W oil provided a wide range of damper 
fluid viscosity. Copper-Constantan thermocouples were used to monitor the 
oil inlet and outlet temperatures from the dampers. A thermal well was used 
for monitoring the pressurized inlet temperature. 

Vibration detection was accomplished with Bently-Nevada non-contacting 
proximity displacement transducers. The output of the proximitors was wired 
into a digital vector filter (DVF2) also manufactured by Bently-Nevada. The 
DVF2 provided a digital readout of peal to peak displacement, RPM, and phase 
angle for location of the rotor's 'high spot' for balancing. Two Tektronix 
oscilloscopes, an X-Y-Y plotter, and an HP spectrum analyzer were also 
useful aids in studying and recording vibraitonal response. The spectrum 
analyzer was expecially useful for investigating nonsynchronous responses. 

A five thermocouple input, digital readout potentiometer was used to monitor 
oil temperature into and out of the SFD. The schematic diagram for this 
assembly is shown in Fig. 33. 

9.3 Experimental Method and Objectives 

The initial purpose of the experimentation was to verify that the test 
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apparatus could work properly as designed and that operational data could be 
taken as desired. Once the system was debugged, the next step was to 
conduct a brief parametric study. 

The parameters that could be contolled for both the 'short' and 'long* 
bearing configurations of the damper were as follows: 

1) rotor speed 

2) unbalance 

3) clearance 

4) oil temperature 

5) oil inlet pressure 

The second purpose was then to discover the effect on the rotor 
system's performance of varying each parameter while holding the others 
constant. Comparisons could then be made with other authors' results and 
with a simplified analysis and conclusions made. 

The third and final purpose was to recommend possible improvements to 
the test equipment and to suggest further work. 

Bode plots of peak-to-peak amplitude and phase vs. rotor speed could be 
made at both damper locations and at any point along the shaft for both x 
and y positions. Also, damper and shaft orbits could be obtained and photo- 
graphed for a particular rotor speed. The orbits (oscilloscope tracings) 
have a gap in them made by a reference mark useful for phase studies. The 
phase mentioned here is the angle between the high spot of the rotor and a 
reference mark on the shaft. The meaningful information is the change in 
this angle with speed (denoting changes of angles between the high spot and 
the unbalance position). 

The first task was to balance the rotor after the system was assembled 
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with the nul 1-clearance inserts in the dampers providing 'rigid' bearing 
supports. The DVF2 aided in this endeavor. Fig. 34 shows the midspan 
vibrational response for the balanced rotor. Note that all Bode plots were 
made from a location two inches from the midspan location. Comparing with 
the other plot on the same Figure for a small unbalance, it can be seen that 
the residual unbalance is small. The sharp critical speed response on Fig. 
34 and the abrupt phase change shown on Fig. 35 for a small unbalance demon- 
strates the high stiffness and lack of damping of the ball bearings. The 
first critical speed of the rotor can be seen to be about 9200 RPM. 

In practice, the four damper mounting screws on each SFD were used to 
adjust the 0-ring squeeze to seal adequately, yet contribute minimal imped- 
ance to journal motion. Shims were then required to restrain the damper 
bearings from vibrating relative to their respective housings. 

When damper orbit studies were made, the GAP button on the DVF2 was 
used to obtain the approximate orbit center coordinates within the clearance 
circle. The GAP button provides a measure of the distance between the 
proximeter probes and the rotor. This method was also used to measure the 
radius of each clearance circle. 

Fig. 36a and b are examples of typical photographs of the damper and 
shaft orbits respectively, taken from oscilloscope tracings. The rotation 
of the shaft is clockwise with respect to all orbits presented in this 
paper. Copies of each orbit photograph of interest were reduced to their 
proper scale, and placed in the correct position and orientation on drawings 
of their particular clearance circle. This was done because knowledge of 
the location (offset from center) of the orbit within the clearance circle 
was found to be useful. 
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9.4 Experimental Remits 

The complete set of experimental results are presented in reference 
[16] for the rotor system supported by SFD without end seals (annular oil 
inlet/outlet) and with end seals (port inlet/outlet). A summary of these 
follows. 

Fig. 37, a superpositioning of four separate damper orbits, shows the 
effects of changing the unbalance while holding the other parameters 
constant. Increasing the unbalance increased the vibration of the rotor 
which resulted in a further 'lifting off of the dampers from the bottom of 
the clearance circle as well as larger elliptical orbits. 

Figs. 38 through 42 show amplitude plots of the rotor's x-axis for 
different values of unbalance. Note that the trends are consistent despite 
the clearance size. All cases show that at least as a* small amount of 
effective damping was present, yielding a response superior to the 'rigid' 
response of Fig. 34. This is apparent from the lower amplitudes and the 
less abrupt phase shifting with additional dissipation. As more unbalance 
was applied, larger orbits were observed and the effective damping 
increased, displaying the amplitude dependent damping characteristics of the 
SFD. As noted by Cookson and Kossa for this type of dampers, the critical 
speed did not change appreciably. This shows that SFD act only as energy 
dissipators not as springs in series with the rotor. Note that the 
resonance peaks are not well defined, appearing as eroded versions of what 
they might have been with 'rigid' supports. The onset of a noise which was 
taken to be that of cavitation was noted at the peaks where the 'erosion' 
begins. Comparing the plots for the various clearances, it can be seen that 
larger clearances yielded freer motion of the damper journals and thus 


'softer' diaper*. Larger clearances yielded more favorable rotor responses 
(for this system). 

5.5 Comparison with ± Simol if ied Computer Model 

The rotor was modelled using an available code, assuming the rotor to 
be rigid [8]. Of coarse, the rotor was not rigid at the speeds which the 
experimental resnlts were obtained. Thus, the damper excitation was not 
necessarily the same and direct comparisons could not be made. However, 
qualitative comparisons concerning the general character of the orbits were 
found to be useful. All computed cases presented were run at 10,000 RPM, 

9 

with 3 gm unbalance and low inlet pressure using a short bearing approxima- 
tion of the Reynolds equation. 

Considering Figs. 43, 44, and 45, with increasing viscosity, the 
orbit's amplitudes and offsets decreased, showing the damper to be growing 
stiffer. This agrees with the experimental results. Also, as the offset 
decreased, the orbits became more circular, which also agrees with 
experimental evidence. All of the computer produced orbits (clockwise ro- 
tations) are left of center (e.g., see Fig. 44) and approach the center line 
(y-axis) with decreasing offset. This is also the general rule for the 
experimentally obtained orbits. From this last comparison it can be seen 
that increasing the clearance in both cases produced larger orbits and off- 
sets and thus a softer damper. 
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10. SUMMARY AND CONCLUSIONS 


Phase -1 of this work provided the following: 

• Overall solution strategy for the general engine dynamics problem 

• Clear definition of where current efforts shonld be concentrated 

• A SFD software package ready for FE implant 
Phase -2 of this work provided the following: 

• Methodology for implanting SFD element into a general FE code 

• Benchmarking the SFD/ADINA combination 

• Numerous relevant example-problem solutions 
Phase -3 of this work provided the following: 

• Conception, design and construction of a SFD/flexible rotor test rig 

• The results of basic experimental parametric studies 

• Comparisons with computed results and results of previous 
investigators 

It is concluded that general engine dynamic analyses as a standard 
design-study computational tool is readily available, practical and highly 
desirable for the prediction and understanding of complex engine dynamic 
behavior. Improved definition of engine dynamic response will likely 
provide valuable information and insights leading to reduced maintenance and 
overhaul costs on existing engine configurations. Furthermore, application 
of advanced engine dynamic simulation methods will provide a considerable 
cost reduction in the development of new engine designs by eliminating some 
of the trial-and-error process done with engine hardware development. 

Clearly, the field of aircraft engine dynamics is presently in a 
position where there is both a need for substantial advances and feasible 
means by which such advances can be accomplished. The necessary approach 
and methods have been developed, proven and demonstrated by this work. 
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NOMENCLATURE 


C *= radial clearance of damper annulus 

D = nominal damper annulus diameter e 2R 

e = damper eccentricity 

F x = X - component of damper force 

F y • I - component of damper force 

h = damper annulus film thickness distribution 

L = damper length 

p *= damper film thickness distribution 

R * nominal damper annulus radius 

t * time 

z « R0 * damper annulus circumferential coordinate 
X ■ X - direction radial motion coordinate 

Y * Y - direction radial motion coordinate 

x * damper annulus axial coordinate 

p ■ damper lubricant viscosity 

C * frequency of vibration excitation 

SPECIAL TERMINOLOGY 

Infinitely Long Bearing Model - axial flow is neglected (-r— <<-|— ) 

oZ o X / 

3 3 

Infinitely Short Bearing Model - circumferential flow is neglected (|— <<|y) 

Driver Code - Any computer code which calls the squeeze-film damper force 
computation code 
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Fig. l Input/Output of Damper Pilot Code 
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Fig. 2 Inertial coordinates 
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Fig. 3 A typical aircraft engine damper configuration 
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Fig. 4(a) Configuration frequently used In military applications 
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Fig. 4(b) Configuration frequently used in coranercial applications 
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Unwrapped squeeze-film pressure solution domains for configurations 
shown in Figure 4. 



(a) e/c*0.05 



(b) e/c-0.20 


Figure 6 Pressure distribution in circumferential direction and time 
of one cycle of circular orbit(long-bearing solution). 


(d) e/c-0.95 
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Pressure distribution in circumferential direction and time 
of 6 one cycle of circular orbit (long-bearing solut 





(a) e/c*0.05 



(b) e/c*0.20 


Figure 7 Pressure distribution in circumferential direction and time 
of one cycle of circular orbit (short-bearing solution). 
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figure 8 


Simple 2-mass ,4-degree 
(Same damper parameters 


freedom. Test case 
as on page |g) 
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(a) Rotor and stator orbits 



(b) Rotor orbit relative to stator 
(clearance circle shown) 


Fig. 9 Nonlinear dynamic transient of simple 4 DOE system(See Fig. 8) 
J F | *1 00 lbs> = 150 rad/sec ,fU=M2=500 1 bs ,Kx=Ky= 1 1 6000 lbs/in. 


52 




(a) Rotor and stato 



(b) 


Rotor orbit rel 
(clearance circ 


Fig. 11 Nonlinear dynamic transient 
| F J =300 1 bs‘/ t ' = 1 50 rad/sec, M 
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M 




— r~ 

!» 


(b) Rotor orbit relative to stator 
(clearance circle shown) 


Fig. 13 Nonlinear dynamic transient of simple 4 DOE system 

|F| *1000 IbsV"* 150 rad/sec ,M1*M2*500 lbs ,Kx=Ky=116000 lbs/in 
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. I5SEC 

SHOCK LOAD 


Fig. 1A. 


S 




Z(IO') 


r\ = | 50 RAD/SEC 
1P| = 20 LBS 
IS 1=0 LB 



cvjj 

» 



£1=150 RAD/SEC 
IPI = 20 LBS 
ISI =0 LB 



Fig. 16. Stator Displacement j 
|P| • 20 lbs. 












ZOO ) 



Fig. 19. Stator Displacement; 
|P| - 2500 lbs. 


02 



•n. = 150 RAD/SEC 
= 2500 LBS 
= 1000 LBS 



-04 


.0 


.0 


FIGURE 6.14 ROTOR 0ISPLACEMEN1 

(SEE FIG. 6.1 ) 


Fig. 20. Rotor Orbit; 

|P| - 2500 lbs. 



A El B El C El D M s =Mj fdb =- 50QLBS 
r*~ L 4 * L -«t*~ L— H M =IOOOLBS 
E = 30.*I0 6 PSI , I =4 63. INCH 4 , L = 48INCH 
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CVJ 




-08 



n=750 RAQ/SEC 
P =2500 LBS 

S =0 LB 


Fig. 24. Rotor Displacement; 
|P| - 2500 lbs. 

0 * 750 rad/sec 
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/\ = 1 50 RAD/SEC 

I P I = 500 LBS 
|S |* O' LB 



r 


Y(IO ) 

.0 .12 


Fig. 25. Rotor Orbit; 

| P | = 2500 Lbs. 

0 =750 rad/sec 





Figure 27. Experimental lest Facility Rotor 


original PAGE IS 
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Figure 28 . Schematic Side View of Experimental Damper 
Bearing Configuration 
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1/2 Scale 
A1 umi num 
2 Reqd. 


Fig. 29. Damper Bearing 






R, * 1-5/32" 
R 2 « 1-19/32 
R 2 - 2-5/8" 


1/2 Scale 
A1 uml num 
2 Reqd. 


Bearing Housing Side View 







Pi gure 34 


PHASE 


Midshaft Phase vs. Rotor Speed 
Rigid Supports 
u 8 .2 gm 



Figure 35 
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b. Shaft Orbit, same as above except 2 psl, 2.5 MIl/DIV. 

Figure 36. Photographs of Oscllllscope Tracings of Damper 
and Shaft Orbits 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Damper Orbits Relative to Clearance Circle 
Groove Feed, 80°F, P in * 2 psi, 9000 RPM 


Y, 



Figure 37 
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Figure 38 
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PHASE 


Midshaft Phase vs. Rotor Speed 


Groove Feed, 80°F, C * .010% P in « 2 psi 

1 U * .0662 (1 gm) 

2 U « .1987 (3 gm) 



Figure 39 
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Midshaft Displacement vs. Rotor Speed 
Groove Feed, C » .010", 80°F, Pj - 2 psi 



Figure 40 
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PHASE 


Hi dsh aft Phase vs Rotor Speed 

Groove Feed, 80°F, C - .006", p « 2 psi 

1 u * 1 gm 

2 u ■ 2 gm 



Figure 41 
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Damper Orbits Relative to Cl earance Circle 

Groove Feed, Low P lRf u * 3.2 gn. 10,000 RPM 

1 y * 10’ 6 Reyns. ( 1 6 0 0 F ) 

2 y * 5x1 0” 6 Reyns. (100°F) 

3 y « 10x1 0" 6 Reyns. (80°F) 



Figure A3 
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Damper Orbits Relative to Clearance Circle 


Groove Feed, Low P 1n , u * 3.2 gm, 10,000 RPM 

1 10"* Reyns. (160°F) 

2 5x10"* Reyns. { 1 00 ° F ) 



Figure 44 
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Damper Orbit Relative to Clearance Circle 


Groove Feed, Low P^ n » u * 3.2 gm, 10,000 RPM 

1 10’* Reyns. ( 1 60 °F ) 

2 10x1 0” * Reyns. (80°F) 


Y 



Figure 45 
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